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Abstract 

Turbulent flow fields can be expanded into a series in a set of basic 
functions. The terms of such series are often called modes. A master- 
(or determining) mode set is a subset of these modes, the time history of 
which uniquely determines the time history of the entire turbulent flow 
provided that this flow is developed. In the present work the existence 
of the master-mode-set is demonstrated numerically for turbulent channel 
flow. The minimal size of a master-mode set and the rate of the process 
of the recovery of the entire flow from the master-mode set history are 
estimated. The velocity fleld corresponding to the minimal master-mode 
set is found to be a good approximation for mean velocity in the entire 
flow field. Mean characteristics involving velocity derivatives deviate in a 
very close vicinity to the wall, while master-mode two-point correlations 
exhibit unrealistic oscillations. This can be improved by using a larger 
than minimal master-mode set. The near- wall streaks are found to be 
contained in the velocity fleld corresponding to the minimal master-mode 
set, and the same is true at least for the large-scale part of the longi- 
tudinal vorticity structure. A database containing the time history of a 
master-mode set is demonstrated to be an efficient tool for investigating 
rare events in turbulent flows. In particular, a travelling-wave-like object 
was identified on the basis of the analysis of the database. Two master- 
mode-set databases of the time history of a turbulent channel flow are 
made available online at http://www.dnsdata.afm.ses.soton.ac.uk/, The 
services provided include the facility for the code uploaded by a user to 
be run on the server with an access to the data. 
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1 Introduction 
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Direct numerical simulations of turbulent flows are now quite common. They 
have proved to be an important tool for fundamental research on turbulence, 
although calculating the entire turbulent flow around an aircraft is not yet 
possible. Indeed, even for moderate Reynolds numbers and simple geometries 
the computational effort required for direct numerical simulation is very high. 
For example, the recent simulation of a plane channel flow for Rcr = 2003 
{Rct is the Reynolds number based on friction velocity) required about 6 x 10^ 
processor- hours ( iHoyas &: Jimenezl . 120061 ) . Very few research groups have access 



to computational resources of that scale but even for them repeating a major 
direct numerical simulation is not always possible. It is natural, therefore, to 
store the data for future analysis, and indeed, databases of direct numerical sim- 
ulation of turbulent flow are also becoming common. A turbulent flow database 
gives a considerable advantage, as the data can be quickly accessed at a mod- 
erate cost. However, a database of the full time history of a turbulent flow 
produced by a major direct numerical simulation would be very large. For ex- 
ample, the above-mentioned simulation gave about 25 TB of raw data. Opening 
access (say, via Internet) to such a large volume is not feasible yet. For this rea- 
son the existing open-access direct numerical simulation databases contain only 
a limited number of samples of the instantaneous flow fleld complemented with 
various averaged data. Because of this the majority of researchers working on 
turbulence have a limited access to full direct numerical simulation results. 

Direct numerical simulations are commonly performed for flows in bounded 
domains (say, a pipe flow periodic in axial direction). It is generally agreed 
that in this case as time advances the state of the flow approaches a flnite- 
dimensional attr actor. While the dimension of this attractor is believed to be 
high for high Reynolds numbers, it can be expected to be much smaller than 
the number of modes required for a fully resolved direct numerical simulation. 
Therefore, at least for developed turbulent flows it should be possible in principle 
to store only part of the modes used in the direct numerical simulation and still 
have a complete representation of the flow. This could resolve the problem 
of an efficient use of direct numerical simulation results by a wide research 
community. At the same time a study in this direction can give the answer to 
a deep fundamental question concerning the nature of turbulence, namely, the 
question about how much redundant information is indeed contained in a full 
time history of a turbulent flow. 

Determining the dimension of the attractor ( iKeefe et al\ . 119921 ) and stor- 



ing the amplitudes of the number of modes slightly larger than this dimension 
does not actually resolve the issue of the size of a direct numerical simulation 



^Experience of disseminating these results shows that for those who are new to the subject 
understanding can be achieved quicker if it is first pointed out that master modes are not 
"something like POD". In fact, POD modes can also be divided into a master mode set and 
remaining modes. It is also worth mentioning that master modes and determining modes are 
the same: we use the term master modes because we also need the complementary term slave 
modes. 
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database. Indeed, since the attractor itself is unknown, even if the stored modes 
determine the state of the dynamical system on the attractor uniquely, there 
is no way to determine the amplitudes of the remaining modes. It is also not 
clear which modes should be stored. A more appropriate solution is to store 
a set of amplitudes of the so-called master modes, also known as determining 
modes. Let {(pni^)} be a full basis so that a Navier-Stokes solution UA{t, x) can 
be expanded as 



ua = ^ An(t)(^„(x). 



Tl=l 



For double-periodic two-dimensionaQ flows iFoias &: Prodil (119671 ) proved that 



there is a finite such that if another solution 

oo 

ub = ^B^(t)(/)^(x) 



n=l 



has the same amplitudes of the first N modes, that is if A^(t) = B^(t) for 
1 < n < A then the two solutions converge with time, that is 



||uyi — u^ll as t ^ oo. (1) 

This implies that for a developed flow knowing only the time history of these 
A amplitudes is sufficient for uniquely determining all the amplitudes. In other 
words, these A master modes provide a full information about the developed 
turbulent flow. The scope of the study described in the present paper includes 
clarifying further this idea as applied to numerical calculations, calculating the 
master-mode set for 3D channel flow, investigating some of the properties of 
this set, developing an online access to a database of the master-mode set, and 
giving an example of its use. 



2 Master- mode set of a numerical solution 
2.1 Definition 

In direct numerical simulation the solution is obtained by marching in time with 
discrete time steps. At each step the solution can be represented as 

s 

Ui(x) = ^Ui,„^^(x), 

n=l 

where the subscript i corresponds to the solution at time ti, and (^n(x) is a given 
set of functions. Depending on the particular method, this set of functions can 

■^For tree-dimensional flows obtaining similar results depends on the ability to prove that 
the vorticity is bounded, which is also needed for proving the existence of the tree-dimensional 
solution of the Navier-Stokes equations (|Constantin et all . [1985 ): this is a famous unsolved 
problem. However, for solutions obtained numerically a further step forward towards rigorous 
results can be made, see (jChernvshenko et g/.l . [20071 ). 
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be a part of a full basis, as it is typical for Galerkin methods, and as is the case 
in the particular numerical method used in the present paper, or just an element 
of some interpolation formula, as in finite-difference methods based on a control 
volume approach. Each term in this sum is a mode, and u^^^ is the amplitude 
of the mode n at the step i. Then one time step can usually be described by the 
formula 

Ui+l,n = ^n(Ui,l, Ui,2, . . . , Ui,s), 72 = 1, . . . , 5", (2) 

where the functions depend on the particular numerical scheme. (Note 
that we consider the case when the boundary conditions and body forces are 
independent of time so that the functions Dn are independent of i. Note also that 
in some schemes Ui_|_i depends not only on but also on u^-i. For simplicity we 
exclude such schemes from our analysis.) Let u^^^, z = 1, 2, . . . be a numerical 
solution satisfying ([2]). For a given selection n = mi,m2, ....rriK of the mode 
numbers consider a sequence {vi(x)} : 

5 

Vi(x) = ^ Vi,„^„(x) 
n=l 

such that 

t+i,n I i}^^(v-i,Vi,2,-..,Vi,5), n^rrij, j = l,...,K. 

In other words, values of Vi_|_i „^ are obtained with the time-marching formula 
([2]) applied to v^^^ if they do not belong to the set of modes in question and 
are taken to be equal to the corresponding modes of the numerical solution 
otherwise. 

If for any Vi(x) 

I |vi(x) — Ui(x) 1 1 ^ as i ^ oo (3) 

then we will call the set M of modes with n = mi, m2, rriK the master-mode 
set of the numerical solution Ui(x). It should be noted that the sequence Vi(x) 
does not approximate a solution to the Navier-Stokes equations, in contrast 
to ([I]) where ua and are both solutions of the Navier-Stokes equations. It 
also should be noted that the selection n = mi, m2, giving a master- 
mode set for one solution does not necessarily give master-mode set for all other 
possible solutions of the same problem but with different initial conditions. 

There is always a trivial master-mode set including all the modes of the 
numerical calculation. In this case K = S. The natural goal, however, is to 
find master-mode sets with as small a, K a,s possible. It may well be that the 
minimal master-mode set is not unique: the property ([3]) is the property of the 
set of modes and not the property of individual modes. For this reason we use 
only the term master-mode set rather than the term master mode in the rest of 
the paper. 
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2.2 Master code and slave code 



For any given set M of modes there is an obvious way of verifying if this set is 
a master- mode set. The method consists in simultaneously running two direct 
numerical simulation codes, which it is natural to call the master code and the 
slave code. The master code is just a standard direct numerical simulation code 
calculating the solution The slave code is an exact copy of the master 

code except that at each time step the amplitudes of the modes belonging to 
M are replaced by their values calculated by the master code. If the initial 
conditions for both codes are identical, the entire master and slave solutions will 
be identical, too. If the initial conditions are different then the slave solution 
will tend to the master solution when set M is a master-mode set. Otherwise, 
the solutions will deviate. 

If the modes are ordered by the likelihood of them belonging to the master- 
mode set then the minimal master-mode set can be found by dichotomy in the 
set size K. This approach was used by lOlson fc TitH (l2003h . iHenshaw et"al] 
( 120031 ) ■ and lYoshida et all ( 120051 ) in their studies of isotropic flows and proved 
to be highly efficient. In general, the result of such a calculation depends on 
the method of ordering the modes. In the studies on isotropic flows Fourier 
expansions were used and the modes were ordered by the magnitude of their 
wavenumbers. For a 3D channel flow a Fourier expansion in the wall-normal 
direction is not convenient, and the question of ordering the modes requires 
more attention, as described in the following section. 



3 Calculations of master-mode sets for 3D tur- 
bulent channel flow 

3.1 Problem formulation 

Numerical calculations were performed for the flow in a plane channel. The 
governing equations of incompressible flow, that is the continuity equation and 
the momentum equations, were made non-dimensional with the channel half- 
width h* and the friction velocity u* — \/t^/p*: where is the mean wall 
shear stress and p* is the density. Here asterisks are used to denote dimensional 
quantities. The Reynolds number is Rcr = u*h*/u*. The non-dimensional 
quantities are then deflned as u = u*/u*, x = x* /h*, p = p*/p*u*'^, t = 
t*u*/h*. Using u* as the velocity scale results in the normalised mean pressure 
gradient dP/dx = — 1. 

The continuity equation has the form 



V-u = 



while the momentum equations are 

9u „ „ dP 1 



^ + u.Vu=-Vp-e.^-f^VV (4) 
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where is the unit vector of the x-axis. In the above equation the pressure 
has been spht into mean P and fluctuating p components. The coordinates are 
set as X in the streamwise direction, y in the spanwise direction and z in the 
waU-normal direction, with the channel waUs at z = ±1. No-shp and imper- 
meabihty boundary conditions are imposed on the waUs. Periodicity conditions 
are used in the directions paraUel to the waUs: u(x + L^, y, z, t) = u(x, y, z, t), 
p{x + Lx,y,z,t) =p{x,y,z,t), u{x,y + Ly,z,t) = u(x,y,z,t), and p{x,y, z,t) = 
p{x,y + Ly, z,t), where and Ly are the periods equal to the dimensions of 
the computation domain. 



3.2 The numerical method 

We used a version of the parallel pseudo-spectral channel flow code of lSandham fc HowardI 



( 119981 ) further modifled for calculating master-modes. The modiflcation con- 



sisted in replacing parts of the code used in iChernyshenko fc Baigj (|2005l ) for 



calculating passive scalars with what effectively is another copy of the core 
Navier-Stokes solver, and organising the mechanism of replacing at each time 
step the amplitudes of a certain set of modes in the second solver by their 
amplitudes taken from the flrst solver. In order to ensure that no errors were 
introduced during the modification both solvers were then subject to all the 
standard tests for codes used for direct numerical simulation of turbulence. 

The ISandham fc HowardI (fT998h code uses Fourier series expansion in the 
wall-parallel directions and Chebyshev series in the wall-normal direction, so 
that the solution is represented in the form 

Ny/2 N,-l 
ka; = -N^/2 ky = -Ny/2 k^=0 

where T^^ (z) is a Chebyshev polynomial of order and Uk^^ky,k^ (t) are complex- 
valued amplitudes of the modes. Since the solution is real-valued, U/j k k — 
U-k^,-ky,k^- Continuity imposes additional constraints on the components of 
the vectors u.k^,k ,k^- In what follows, when the number of modes is given we 
will refer to the number of real-valued vectors: for example, with this way of 
counting modes the total number of modes in ([5]) is A^^ x A^^ x Nz- The table 
gives the parameters of the four main calculations referred to in the text, with 
some results described in more detail in what follows. 
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3.3 Mode ordering 

Finding the smallest master-mode set by trial-and-error is inefficient. If modes 
are ordered, a more efficient dichotomy algorithm can be used, or simply a set 
of possible sizes of the master-mode set can be tried. Ordering the modes by 
the squared magnitude k"^ + ky + k^ of their wavenumber is non-physical, since 
this method does not take into account the anisotropy of the flow and since k^, 
being the order of the Chebyshev polynomial, should not actually be directly 
compared with k^ and ky. However, since this method is very straightforward 
it was tried. 

A more physical approach could be based on the averaged energy content 
of the modes. However, the modes (E]) are not orthogonal with respect to the 
standard energy norm. Instead, we introduce a Chebyshev-weighted energy 



norm: 




u 



|u||^ = / // —^===dxdydz. 



2 

~2 



This norm of the solution is equal to the sum of the norms of the modes, and 
hence the modes can be ordered by the time-averaged values of their amplitudes 
{\'^k^,ky,k^ P)- The calculations were run for a certain period of time and statistics 
for calculating these averages was collected. Much longer calculation would be 
required for obtaining accurate values of the averages, but for our purpose high 
accuracy is not required. The modes were then ordered by (|u/c^,fc ,A;^P) with 
the assumption that more energetic modes are more likely to belong to the 
master- mode set. 



Schmidtmann et al\ (119971 . 119981 ) in their work on non-linear Galerkin meth- 
ods suggested that ordering the modes by their enstrophy content may be more 
efficient, and ordering the modes by their weighted enstrophy was also attempted 
in the present study. 

Calculations performed for Rcr = 180 in the box of the dimensions Lj; x 
LyXLz = 4x3x2 showed that the minimal size Nmin of the master-mode set 
exceeds 5000 if the modes are ordered by the wavenumber, that Nmin ~ 4800 
for enstrophy-based ordering, and that 2800 < Nmin < 3500 for energy-based 
ordering. The energy-based ordering was then used in all other calculations. 



4 Properties of the master-mode set 
4.1 Minimal size of a master- mode set 

The most widely-known estimate of the number of degrees of freedom of a tur- 
bulent three-dimensional flow follows from the arguments of iLandau &: Lifshitzl 
(fl959D : 



N ~ {Lo/Ldf, (6) 
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64 X 64 X 64, Nmm = 1600, 2800, 2900, 6000, 19000. 



where Lq is the typical large lengthscale and is the Kolmogorov dissipation 
length given by 

La = v^l'le'l\ 

Here, v is the kinematic viscosity and e is the energy dissipation rate per unit 
Constantin ei aL\ (1198.^ gave a rigorous mathematical derivation of ([6]) 



mass. 

as an upper bound of the number of degrees of freedom. They also pointed 
out that in order to deduce the dependence of on the Reynolds number, an 
assumption is needed about the shape of the energy spectrum. Assuming the 
Kolmogorov law of 5/3 leads to the widely accepted estimate 

N ~ Re^/^. (7) 

However, in general the spectrum can have such a shape that the estimate 
becomes 

N ~ Re^. (8) 

They then point out that since the Kolmogorov law is obtained for homogeneous 
isotropic high-Reynolds-number turbulence, ([7]) can be invalid for relatively- 
low- i?e, inhomogeneous flows characteristic for direct numerical simulations of 
turbulence, and express a concern about the validity of numerical simulation for 
which the number of degrees of freedom is often predicted on the i?e^^^ rule. 
On the other hand, iHenshaw et al\ (120031 ) argue that the number of degrees 



of freedom can be even less than ([7]). They considered the case of decaying 
isotropic flow, which is difficult to compare directly to the developed channel 
flow considered in the present study. Their main idea (expressed using our 
terminology) is that a master-mode flow fleld can form regions of high dissipation 
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Figure 2: Convergence history: Rcr = 360, Lj^ = 6, Ly = 3, AT^. 
128 X 128 X 160, Nmm = 18365,22879,29716,34601,41069,63448, 



X Ny X 



(shear layers or concentrated vortices) and it can advect the errors in the slave- 
mode field to these high-dissipation regions where the errors will quickly decay. 
Both the advection rate and the decay rate may be even independent of Re (the 
latter because as Re increases the thickness of, say, shear layers can decrease 
accordingly). This was demonstrated to be sometimes the case for Burgers 
equations. For the Navier-Stokes equations, however, iHenshaw et al\ ( 120031 ) 



observed an increase in the number of degrees of freedom as viscosity decreased, 
but the rate of this increase was not evaluated. 



Yoshida et all (120051 ) considered master modes in a forced isotropic flows and 



their calculations confirm d?]). However, this does not refute the lConstantin et al. 



(119851 ) argument, since one of the possible reasons for the deviation from ([7]) is 
the flow anisotropy. 

Calculating the minimum size of the master-mode set is an efficient way of 
estimating the number of the degrees of freedom in the turbulent fiow, thus 
addressing the question of which rule, (I7j) or ([8]), is the more appropriate for 
the channel fiow. Figures [1] and l2] show the time-dependence of the L2 norm of 
the difference between the master-code solution u.^^ and the slave-code solution 
Us for different sizes A^mm of the master-mode set for two of the calculated 
cases. The calculations were performed with single precision, so that the random 
oscillations at the bottom of figure [T] are the result of the rounding error noise. 
In figure l2] the effect of rounding errors is averaged over a larger volume and 
oscillates less. These figures give the bounds for the minimum size of the master- 
mode set, namely, for Rer = 180 case 2800 < Nmin < 2900 (figure [1]), and for 
Rer = 360 case 22879 < N^i^ < 29716 (figure E]). 

Note that the convergence is not monotone for smaller Nmm and, hence, the 
convergence rate (determined from the slope of the dashed lines in the following 
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analysis) can be found only approximately. 

Figures [T] and [2] establish the existence of the master-mode sets in three- 
dimensional turbulent flow in a plane channel. 

We also checked that increasing the number of modes used in the numerical 
calculations above the necessary minimum does not affect the size of the minimal 
master-mode set. Namely, calculations for Rcr = 360, = Q, Ly = 3 were 
performed with an increased number of modes 256 x 256 x 160. The minimal 
size of the master-mode set was found to be about 29000, which is the same as 
the minimal size of the master-mode set for 128 x 128 x 160 modes calculation. 

Two other calculations were also performed, for Rcr = 85, Lj^ = 5, Ly = 5, 
iV^ X A^y X A^^ = 32 X 32 X 32, which gave 595 < Nmin < 937, and for Rcr = 360, 

= 14, Ly = 9, X Ny X = 512 x 512 x 160, which gave 161336 < 
Nmin < 225900. The calculations with Rer — 85 were performed for the sake 
of comparisons with the only other known previous attempt to calculate the 
number of the degrees of freedom of a plane channel flow. iKeefe et al\ ( 119921 ) 
obtained an approximate estimate Dx ^ 780 of the Lyapunov dimension Dx 
of the attractor for a plain channel flow in a box of the size Lj; = Ly = 1.67r 
using two calculations with two grids, 16 x 33 x 8 and 32 x 32 x 32. For the 
latter grid only partial results were obtained. We had to use a slightly higher 
Reynolds number since at Re^ = 80 we were unable to obtain a turbulent flow. 
Note that the number of degrees of freedom corresponds to twice the size of 
the minimal master-mode set since in our notation each mode represents two 
degrees of freedom, the mode amplitude being a vector the three components 
of which are related by the continuity equation. 

The Lyapunov dimension of the attractor and the minimal size of the master- 
mode set are not the same but they both provide a measure of the number of 
the degrees of freedom of the flow. In view of this and all other factors the 
agreement between IKeefe et al.\ (119921 ) and our results is quite good. The case 
for Rcr = 360, Lj; = 14, Ly = 9 was selected as the largest calculation for which 
we could produce a master-mode set time-history database with the available 
computational resources. 

Obviously, the number of the degrees of freedom of the flow is dependent 
on the dimensions of the computational box. The easiest and most natural 
assumption is that a factor L^Ly should be included in the right-hand-side of 
(17|) and dHj). Note that L^; and Ly are the non-dimensional ratios of the length 
and depth of the computational box to the half-width of the channel. Also, 
Re in these formulae should be based on the maximum velocity (see (28) in 
Constantin et all ( 119851 )). Hence, flgure I3] shows a plot of Nmin/iLxLy) versus 



Re based on the maximum velocity for all available data together with lines 
corresponding to d?]) and dH]). (The IKeefe et al\ ()1992l ) Re-r — 80 point is on 
the right from our Rer = 85 point because the maximum mean velocity in 
their calculation was somewhat greater than in our calculations.) The avail- 
able data flt (IE]) for smaller Re and (I7j) for larger Re. This suggests that the 
Constantin et al\ (119851 )) argument may be correct for smaller i?e, for which it 



was proposed, while the Kolmogorov spectrum is emerging at larger Re. More 
data is needed for full certainty. 
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Figure 3: Error bounds on the size of the minimum master-mode set for cal- 
culated cases and two estimates, Re^^^ (solid line) and Re^ (dashed line). 
Rcr = 85, L-c X = 5 X 5, o; Re^ = 180, L^; x = 4 x 3, □ (error bounds too 
close to be seen); Rcr = 360, L^; x = 6 x 3, A; Rcr — 360, L^; x Ly = 14 x 9, 



(overlaps with the previous case); iKeefe et all (119921 ) 



4.2 Convergence rate 

The rate at which the difference between the master solution and the slave 
solution tends to zero can be estimated in the following way. Let Nmm be the 
number of the modes whose amplitudes are being passed from the master to 
the slave code. When Nmm > Nmim the difference between the master and the 
slave solution converges exponentially, Hu^ — Us|| ~ e~"* ( 1 Olson &: Titil . [20031 ) . 



The value of a depends on Nmm- Assuming that the decay is due to viscous 
effects gives an estimate a ~ L^Re~^ , where Lm is the largest characteristic 
length scale of the slave mode set (non-dimensional quantities corresponding to 
(jl]) are used throughout). One can assume that Lm is proportional to the cubic 
root of the volume of the computational box divided by the number of modes: 
Lm ~ coivst{LxLyLz/NmmY^^ ■ This gives the formula 

C ( Nmm ^ ^ 



R^T V-^x-^y-^zj ^ 

It is equivalent to a*L*m ji' = a* {Lx* Ly* L*/Nmm)'^^^ = const, where * de- 
notes dimensional quantities and u is the kinematic viscosity. Naturally, as Nmm 
approaches Nmim ol should tend to zero. This is not reflected by (I9j), which, 
therefore, can be expected to describe an asymptote approached from below as 

Nmm OO. 

Figure m shows the values of aRtr plotted against Nmm/ [LxLyLz] for all 
calculated cases together with the line given by (ED with C = 30. The values 
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Figure 4: Convergence rate, L^ = 2 : Rcr = 85, L^; x L„ = 5 x 5, o; Re^ = 180, 
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of a were obtained by estimating the slope of the plot of the L2 norm (that 
is the square root of the integral of \um — u^p over the calculation domain) 
versus time. It appears that 1^ can indeed be used for rough estimates of the 
convergence rate at large Nmm- 

Yoshida et al\ ( 120051 ) found that the convergence rate multiplied by the Kol- 



mogorov time scale and plotted versus the cut-off wavenumber multiplied by 
the Kolmogorov length is independent of the Reynolds number. In a chan- 
nel flow the energy dissipation rate depends on the distance to the wall. We 
can, however, introduce an average energy dissipation rate (per unit volume) 
as e = Umeam where Umean IS the mean velocity (note that due to the scaling 
used in our work the pressure gradient is equal to —1, and all the quantities are 
non-dimensional). We then introduce the Kolmogorov time and length scales 
as r = {eRer)~^^'^ and rj = (ei?e^)~^/'^. In lYoshida et al\ ( 120051 ) the number of 
modes per unit volume corresponding to the cut-off wave number ka can be eval- 
uated as (4/3)7r/c^/(27r)"^ while in our case it is Nmm/{LxLyLz). Accordingly, 
we deflne our cut-off wavenumber as 



kn = 2n 



3N„ 



4tvL, 



LyL^ 



1/3 



Figure [5] shows the dependence oi a ~ ar on kaf]. The dash-dotted line is an 
approximation of the lYoshida et al\ (120051 ) results. Of course this comparison 
is only tentative since the channel flow is not isotropic and homogeneous. How- 
ever, our results at least do not contradict the conclusion made in lYoshida et al. 



(120051 ) that the Kolmogorov-scaled convergence rate is determined by the quan- 
tities for the energy dissipation scale rather than the energy containing scale. 
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Figure 5: Convergence rate. The dash-dotted hne is an approximation of the 



Yoshida et al\ (120051 ) results. Other hues and symbols are the same as in figure IH 



In large eddy simulation approach the full solution is also divided into large 
scale and small scale components. Then the large scale motion is simulated nu- 
merically with the effect of small scale on large scale described with a subgrid- 
scale model. If a similarity between the large scale motion in the large eddy 
simulation approach and master-modes is assumed then from the results for mas- 
ter modes it follows that subgrid scale models should have memory, for example 
in the form of relaxation relationships. The above data on the convergence rate 
can be used to estimate the time constants of these relaxation relationships. 



4.3 Mean flow 

A master-mode set contains full information about the particular realisation of 
the developed turbulent fiow, since there is a method for recovery of the rest of 
the modes of that realisation. This makes the master-mode set a natural choice 
for storing the time history of the flow in a database. Of course, the velocity 
field corresponding to the master-mode set alone does not exactly coincide with 
the flow. However, for applications where the approximate representation given 
by the master-mode set itself is sufficient, using it is convenient and gives a 
considerable saving of computing resources. Figure [6] shows the comparison 
of the mean velocity profile obtained from the master-mode set alone and the 
velocity profile of the full (that is obtained from all the modes used in the direct 
numerical simulation) fiow for Re = 360, = Ly = 3 and the master- mode 
set size A^™™ = 29716, which is close to the minimal master-mode set size. 
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Figure 6: Mean velocity, Re^- — 360, L^: x = 6 x 3, full flow (solid curve) and 
a minimal master-mode set (Nmm = 29716), A. 

4.4 Turbulent kinetic energy budget 

For a channel flow 



notation is used with x = xi, y = X2, z = x^, u = ui, v = U2, and w — u^, with 
summation over repeated indices. 

Figure [7] shows the various components of the turbulent energy budget for 
the same calculation and the same master-mode set. The pressure term 11 is 
not shown since it is small, however for 11 the agreement between the minimal 
master-mode set field and the full flow is very good. Note the use of the log- 
arithmic scale. There is a significant deviation of the master-mode set curves 
from the full flow results. This deviation, however, is limited to a few near- 
wall points. The reason for this discrepancy becomes evident from the plot 
of the correlation between the minimal master-mode set velocity and the 
slave- mode set velocity u^,^ (not to be confused with and u^), also given 
in the same figure for the wall-normal velocity component. Away from the wall 
the correlation is small, but near the wall the correlation coefficient is close to 
— 1. Regression analysis shows that close to the wall Usm ~ —'^mm, while the 
full velocity u = -|- is significantly smaller in magnitude than both 

Ujnrn and Usm- This happens because the Chebyshev polynomials, which con- 




= P + T + D + U + ^, 



where P 



and $ 
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Figure 7: Turbulent kinetic energy budget components: full flow (solid curves) 
and master-mode set only (dashed curves); production P, -f; transport T, •; 
viscous diffusion D, ■; viscous dissipation, E A; total budget for the master- 
mode set (thick solid line). Correlation between the master-mode set and slave- 
mode set wall- normal velocity fields (thick dashed curve with x). Rer = 360, 
X = 6 X 3, N^^ = 29716. 

situte the modes, strongly oscillate near the wall, while the full flow velocity 
does not. This is another indication of the well-known fact that Chebyshev 
polynomials over-resolve turbulent flow near the wall if they fully resolve it 
away from the wall. One can obtain the same conclusion using, for example, 
the so-called directional dissipation scale (iManhartl . 120001 ). The amplitude of 
the modes containing lower-order Chebyshev polynomials is determined by the 
solution behaviour over the channel-width scale. The near-wall over-resolving 
property results in the level of derivatives corresponding to these modes being 
well above the magnitude of the derivatives of the solution. Near the wall the 
higher-order modes, therefore, have to cancel out the lower-order modes. This 
phenomenon manifests itself when the solution is divided into a master mode set 
and a slave mode set. Probably, using the high correlation between master and 
slave modes near the wall the agreement can be improved; away from the wall 
the minimal master-mode set provides a reasonable approximation for energy 
budget components. 
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Figure 8: Streamwise {Ay = 0) two-point correlations, curves without symbols, 
full flow; curves with symbols, master-mode-set; longitudinal velocity correla- 
tion. Rum spanwise velocity correlation, Ryy, ▲; wall- normal velocity cor- 
relation, Rww, •■ For Rcr = 360, = 5.6, L^; x = 6 x 3, and minimal 
master-mode set {Nm = 2900.) 

4.5 Two-point correlations 

For any given wall-parallel plane with coordinate 2; in a plane channel two-point 
correlation coefficients can be defined as 



This definition exploits the homogeneity in the wall-parallel x and y directions. 

Figure [8] shows the streamwise [Ay = 0) two-point correlation coefficients 
for the full ffow and the minimal master- mode set velocity field. The fuU-fiow 
and the master-mode-set curves for the longitudinal velocity and the spanwise 
velocity correlations are in a reasonable agreement. However, the correlation 
for the master-mode- set wall-normal velocity exhibits oscillations around the 
respective fuU-fiow correlation for Ax greater than about 75. This can be caused 
by the sharp cut-off in the energy spectrum of the master-mode-set field at 
approximately the same length scal^. The cut-off scale is indicated by a vertical 
dashed line in figure [8] and in figure [H which shows the corresponding energy 
spectrum. The relation between the length scale and the wave number, = 
360Lj;/kj;, follows from ([5]). 

For the spanwise correlation (figure [TOj) the master-mode-set correlation be- 
gins to oscillate around the fuU-fiow correlation at smaller scale, which may 

^This explanation was suggested by Prof. I.P.Castro 
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Figure 9: Streamwise energy 

the curve; master-mode-set, +. For Rcr = 360, 
the minimal master-mode set {Nmm = 29716). 
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Figure 10: Spanwise (Ax = 0) two-point correlations, curves without symbols, 
full flow; curves with symbols, master-mode-set; longitudinal velocity correla- 
tion. Rum B; wall-normal velocity correlation, Rww, •• For Rct = 360, = 5.6, 
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6x3, and minimal master-mode set (A^„ 



29716). 
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indicate that the cut-off level in the spanwise spectrum is at this smaller scale. 
The correlation for the spanwise velocity component is not shown because it is 
very close to the correlation for the longitudinal velocity component. Note that 
the oscillations are observed for the wall-normal component only. This may be 
due to the fact that at this distance to the wall {z~^ = 5.6) the wall-normal 
velocity component is much smaller than the other velocity components: the 
slave modes, being selected by an energy criterion (see section 13.31) may give a 
relatively higher contribution to the wall-normal velocity component. Naturally, 
increasing the number of modes in the master-mode set reduces the oscillations, 
see section m 

4.6 Structures 

While the master-mode set contains full information about the developed tur- 
bulent flow in the sense that its time history uniquely identifies the particular 
realisation of the turbulent flow, it is far more difficult to establish whether the 
master-mode-set velocity field contains all the essential flow dynamics. If this is 
believed to be true then considering which organised structures are observed in 
the minimal master-mode-set and which are observed in the remaining modes 
may help to establish a relative importance of different organised structures. On 
the other hand, if the relative importance of different structures is considered 
to be known, the analysis of the structures in the minimal master-mode- set and 
the remaining modes can be used to estimate the dynamical significance of the 
minimal master-mode set. From the results obtained so far it appears, however, 
that the minimal master-mode set contains almost all the structures. 

Figure [TT] shows the visualisation of the near- wall streaks in a wall-parallel 
(xy) plane = 5.6 and the visualisation of the streamwise vorticity in the wall- 
normal main-flow-parallel (xz) plane. It is clearly seen that all the streaks are 
contained in the minimal master- mode set. The same is definitely true for the 
large-scale vorticity structures, while some of the smaller vorticity structures, as 
indicated by the dashed line boxes, are actually found in the slave-set but not in 
the master-mode set. It would be interesting to follow the time history of these 
small structures in order to establish their dynamical significance; however, this 
is beyond the scope of the present study. The comparisons made so far show 
only that the minimal master-mode-set captures most of the structures. 

5 Master- mode database and an example of its 
use 

5.1 Online database 

The properties of the master-mode-set, namely, its relatively small size and the 
possibility of recovering all other modes if necessary, make the master-mode-set 
suitable for creating a time-history database of a turbulent flow. For such a 
database a compromise should be made between the number of modes recorded 



18 
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d) e) f) 

Figure 11: Streaks shown by the cut-off of the longitudinal velocity at = 5.6, 
a, b, c, and streamwise vorticity in the wall-normal main-flow-parallel (xz) 
plane, d,e,f. Full flow, a and d; master-mode-set, b and e; remaining modes, 
c and /. For Rcr = 360, x Ly = 6x3, and minimal master- mode set 
(iV^ = 29716). 



and the recovery rate a (see ([9])). Increasing the number of modes above the 
minimal master- mode set size also improves the accuracy with which the master- 
mode set approximates the flow. In the course of the present study two master- 
mode-set databases were created, both for the channel flow at Rcr = 360. The 
first database contains a time history of the master-mode-set with = 51508 
modes, which is about twice the size of the minimal master-mode set, for the box 
size Lj^xLy =6x3 and the duration T = 40. The convergence rate given by ([9]) 
for this database is ct ~ 10. Therefore, achieving 1% accuracy of the recovered 
solution requires the slave mode calculation time AT = In 0.01 ~ 0.4. The 
total size of the database is 112.7Gb. A fully resolved database would require 
about 2.6Tb of storage, other things being equal. The second database is for the 
box size L^; x = 14 x 9 with = 2170044, which is about five times the size 
of the minimal master-mode- set, T = 15, and a ~ 160, which gives AT ^ 0.03. 
It uses 1.78Tb of disk storage. A fully resolved database would require about 
15.2Tb of storage, other things being equal. 

As an alternative to recovering the slave modes one could simply store a set 
of instantaneous fiow fields of the original calculations with the time interval AT, 
and recover the time history by running the full Navier-Stokes calculation using 
the nearest stored fiow field as the initial condition. The calculations needed 
for recovering the slave modes are of a similar volume to full flow calculations, 
so, the CPU time required would be of the same order of magnitude in both 
approaches. For the second approach one would need to store T/ AT snapshots, 
which would require only about 3.2Gb storage for the flrst database and 270Gb 
storage for the second database. However, in order to reproduce exactly the 
same flow by running a full Navier-Stokes calculation starting from the stored 
snapshot one would need to use a computer of the same architecture as that on 
which the original calculation was performed, as otherwise rounding errors would 
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Figure 12: Spanwise (Ax = 0) two-point correlations: curves, full flow as calcu- 
lated by lHu et al\ (120061 ). x Ly = 12 x 6; symbols, master-mode-set database 
with Nmm — 2170044 and x Ly = 14 x 9; longitudinal velocity correlation. 
Run, spanwise velocity correlation, Ryy, A; wall-normal velocity correlation, 
Rww, •■ For Rcr = 360 and = 10. 

lead to the deviation of the original and reproduced solution due to its instability. 
In contrast, recovery of the slave modes is a stable procedure. The master- 
mode database has the advantage of the possibility to use it as an approximate 
time history description of the flow without doing any recovery at all. The 
increase in size of the databases as compared to the minimal master-mode set 
size ensures a better approximation of the full velocity field. Figure [12] shows 
the comparison for the spanwise two-point correlations calculated by lHu et al. 
( 20061 ) and those obtained from the larger master- mode-set database. Other 
quantities give at least as good an agreement. Moreover, the increased number 
of modes in the larger database removed the near-wall discrepancy illustrated 
by figure [71 Overall, the larger database gives quite a reasonable approximation 
to the full fiow. 

The relatively small size of the master-mode databases allows opening a free 



online access. The databases are currently at http:/ /www.dnsdata.afm.ses.soton.ac.uk^ 




The system allows the user to upload a Fortran 90 code, compile and run it on 
the server with access to the data through a subroutine library. The user can 
then download the output of the code. This provides additional savings of the 
network traffic and appears to be more convenient than downloading the entire 
database to the remote computer. 
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Figure 13: Travelling-wave-like object shown by the longitudinal velocity fluc- 
tuation at = 6.9. Rcr = 360, L^; x L„ = 6 x 3, Nmm = 51508. 



5.2 Travelling waves 

Using a master-mode-set database is most advantageous when an analysis of the 
same data needs to be repeated many times, as when searching for rare events 
and repeatedly adjusting the search criterion. Looking for travelling waves in 
a developed turbulent flow gives an example of such a procedure. Travelling 
wave solutions attracted recently much attention. For plane channel flow they 
were found by lltano fc Tohl fl25Mh and iWaleffd (l200ll . l2003h . while the most 



recent attempt to flnd travelling waves in the velocity fields produced by direct 
numerical simulation was performed by iKerswell fc Tuttyl ( 120071 ) , where further 



references can be found. Note that as the Reynolds number and the size of 
the computational domain increase the likelihood of finding a travelling wave 
occupying the entire computational domain decreases. Therefore, one has to 
look for travelling waves in sub- volumes of the domain. 

A travelling wave solution has the form u{x, y, t) = \J{x — ct, y, z), where 
c is the propagation speed of the travelling wave. For such solutions 

du ^ du ^ 

dt dx 

Accordingly, one can assume that a volume W is occupied by a travelling wave 
if 



mm 

c 



J J Jw 



w 



dt dx 



dx dy dz 



is small. To make this idea practical one needs, however, to define exactly 
what is meant here by small, that is to find a proper scale for /. Also, the 
value of c minimising this integral depends on W. The procedure for looking 
for a travelling wave consists then in specifying a set of volumes W and the 
scaling for /, searching over a database for a minimum of the renormalized /, 



21 



looking at the outcome, refining intuitively both the volumes and the scaling 
and repeating the process again. This analysis was performed using the smaller 
database. Eventually, the scale for / was chosen to be 



R 




w 



dvi ^ 



dt 



dx dy dz, 



and the entire calculation domain X Ly X = 6 X 3 X 2 was divided into 
rectangular boxes of the size Ax x Ay x Az = 1.5 x 0.15 x 0.1. The minimum 
of I/R was attained at t = 13.75 in the volume < x < 1.5, 2.4 < y < 
2.55, —1 < z < —0.9. Figure [13] shows the visualisation of the longitudinal 
velocity at several consecutive times, with the area in question marked by a 
dashed rectangle. Note that the visual propagation speed in this figure is about 
10, which is greater than the mean flow velocity at this distance to the wall 
(somewhat less than 7) but which is smaller than the velocity c ~ 14.8 at which 
the minimum of I/R was attained (the dashed rectangle moves at velocity c 
in flgure [T3l). By the visual propagation speed we understand the propagation 
speed of the visually identifled point corresponding to the top of the dark-shaded 
arc in the middle of the dashed rectangle in the flgures for t = 13.7 — 13.875. 
It is interesting to note that these velocities expressed in wall units (about 10 
and 15) are respectively close to the propagation speed of velocity and vorticity 
fluctuations and the propagation speed of pressure fluctuations as identifled 
by iKim &: HussamI (119931 ). Determining whether the object identifled here is 



in fact related to travelling waves or whether it is a certain combination of 
generic near-wall structures is outside the scope of the present paper. The 
observations made here conflrm the efficiency of the master-mode-set database 
as a tool for finding rare events in turbulent fiow. The easy online access to 
the database gives an opportunity for any researcher to further investigate the 
object described in this section. 



6 Conclusions 

Master-modes have a wide scope and have not been previously studied for a 
three-dimensional turbulence in a channel fiow. For this reason each of the 
results described in the present paper may be extended and deepened. One can 
now only summarise the outcomes of the present work: 

• Master-mode sets exist for three-dimensional turbulent flow. 

• The minimal size of a master-mode set appears to be proportional to the 
product of the volume of the computational domain and the 9/4 power 
of the Reynolds number Re (based on the maximal velocity) for large Re. 
A quantitative estimate may be inferred from flgure O In a rather ap- 
proximate way the size of the minimal master-mode-set can be estimated 
as about 2% of the minimal number of modes needed to resolve the flow. 
These estimates may be dependent on the basis used in the calculations; in 
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the present calculations it was Chebyshev polynomials in the wall-normal 
direction and Fourier expansions in the wall-parallel directions. 

• For channel flow the rate of convergence at which the L2 norm of the 
difference between the master solution and the slave solution tends to 
zero can be estimated as ?>Q{N /V)'^/^ / Rer for large A^, where is the 
master-mode set size (that is half of the number of independent real- valued 
parameters needed to identify all the mode amplitudes in the master- 
mode set), V is the volume of the computational domain, and Re-r is the 
Reynolds number based on the friction velocity and the channel half- width. 

• The velocity fleld corresponding to the minimal master- mode set is a good 
approximation for mean velocity and some other mean characteristics. 
However, the minimal master-mode set based on Chebyshev polynomials 
in the wall-normal direction produces large deviations of turbulence kinetic 
energy budget components from their exact values in the close vicinity of 
the wall where Chebyshev polynomials over-resolve the flow. Also, there 
are signiflcant deviations for some of the two-point correlations. Therefore, 
it is reasonable to store in a database a master-mode set of a somewhat 
greater size than minimal. 

• When the flow fleld is represented as a sum of a master-mode set and 
the rest of the modes, the near-wall streaks turn out to be contained in 
the master-mode part. The results for the distribution of the longitudi- 
nal vorticity are less certain, with some small-scale structures sometimes 
observed in the rest of the modes. 

• A database containing the time history of a master-mode set is found to 
be an efficient tool for investigating rare events. In particular, a travelling- 
wave-like object was identified on the basis of the analysis of the database. 

• Two master-mode-set databases of the time history of a turbulent chan- 
nel fiow are available online at http:/ /www. dnsdata.afm.ses.soton.ac.uk/. 
The code uploaded by a user can be run on the server with an access to 
the data. 
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